CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development.

Loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. Here, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. By engineering human induced pluripotent stem cells (hiPSCs) consisting of wild-type (WT), heterozygous (CDX2-Het), and homozygous null CDX2 (CDX2-KO) genotypes, differentiating these cells in a 2D gastruloid model, and subjecting these cells to single-nucleus RNA and ATAC sequencing, we identify several pathways that are dose-dependently regulated by CDX2 including VEGF and non-canonical WNT. snATAC-seq reveals that CDX2-Het cells retain a WT-like chromatin accessibility profile, suggesting accessibility alone is not sufficient to drive this variability in gene expression. Because the loss of CDX2 or TBXT phenocopy one another in vivo, we compared differentially expressed genes in our CDX2-KO to those from TBXT-KO hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and reveal pathways that may underlie the defects in vascular development and allantoic elongation seen in vivo.

TBXT dose sensitivity and the decoupling of nascent mesoderm specification from EMT progression in 2D human gastruloids.

In the nascent mesoderm, TBXT expression must be precisely regulated to ensure that cells exit the primitive streak and pattern the anterior-posterior axis, but how varying dosage informs morphogenesis is not well understood. In this study, we define the transcriptional consequences of TBXT dosage reduction during early human gastrulation using human induced pluripotent stem cell models of gastrulation and mesoderm differentiation. Multi-omic single-nucleus RNA and single-nucleus ATAC sequencing of 2D gastruloids comprising wild-type, TBXT heterozygous or TBXT null human induced pluripotent stem cells reveal that varying TBXT dosage does not compromise the ability of a cell to differentiate into nascent mesoderm, but instead directly influences the temporal progression of the epithelial-to-mesenchymal transition with wild type transitioning first, followed by TBXT heterozygous and then TBXT null. By differentiating cells into nascent mesoderm in a monolayer format, we further illustrate that TBXT dosage directly impacts the persistence of junctional proteins and cell-cell adhesions. These results demonstrate that epithelial-to-mesenchymal transition progression can be decoupled from the acquisition of mesodermal identity in the early gastrula and shed light on the mechanisms underlying human embryogenesis.

CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development.

Proper regulation of gene dosage is critical for the development of the early embryo and the extraembryonic tissues that support it. Specifically, loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. In this study, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. We generate an allelic series for CDX2 in human induced pluripotent stem cells (hiPSCs) consisting of WT, heterozygous, and homozygous null CDX2 genotypes, differentiate these cells in a 2D gastruloid model, and subject these cells to multiomic single nucleus RNA and ATAC sequencing. We identify several genes that CDX2 dose-dependently regulate cytoskeletal integrity and adhesiveness in the extraembryonic mesoderm population, including regulators of the VEGF, canonical WNT, and non-canonical WNT signaling pathways. Despite these dose-dependent gene expression patterns, snATAC-seq reveals that heterozygous CDX2 expression is capable of inducing a WT-like chromatin accessibility profile, suggesting accessibility is not sufficient to drive gene expression when the CDX2 dosage is reduced. Finally, because the loss of CDX2 or TBXT phenocopy one another in vivo, we compare differentially expressed genes in our CDX2 knock-out model to those from TBXT knock-out hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability, including ANK3 and ANGPT1. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and suggest these genes may underlie the defects in vascular development and allantoic elongation seen in the absence or reduction of CDX2 in vivo.

A disrupted compartment boundary underlies abnormal cardiac patterning and congenital heart defects.

Failure of septation of the interventricular septum (IVS) is the most common congenital heart defect (CHD), but mechanisms for patterning the IVS are largely unknown. We show that a Tbx5+/Mef2cAHF+ progenitor lineage forms a compartment boundary bisecting the IVS. This coordinated population originates at a first- and second heart field interface, subsequently forming a morphogenetic nexus. Ablation of Tbx5+/Mef2cAHF+ progenitors cause IVS disorganization, right ventricular hypoplasia and mixing of IVS lineages. Reduced dosage of the CHD transcription factor TBX5 disrupts boundary position and integrity, resulting in ventricular septation defects (VSDs) and patterning defects, including Slit2 and Ntn1 misexpression. Reducing NTN1 dosage partly rescues cardiac defects in Tbx5 mutant embryos. Loss of Slit2 or Ntn1 causes VSDs and perturbed septal lineage distributions. Thus, we identify essential cues that direct progenitors to pattern a compartment boundary for proper cardiac septation, revealing new mechanisms for cardiac birth defects.

TBXT dose sensitivity and the decoupling of nascent mesoderm specification from EMT progression in 2D human gastruloids.

In the nascent mesoderm, levels of Brachyury (TBXT) expression must be precisely regulated to ensure cells exit the primitive streak and pattern the anterior-posterior axis, but how this varying dosage informs morphogenesis is not well understood. In this study, we define the transcriptional consequences of TBXT dose reduction during early human gastrulation using human induced pluripotent stem cell (hiPSC)-based models of gastrulation and mesoderm differentiation. Multiomic single-nucleus RNA and single-nucleus ATAC sequencing of 2D gastruloids comprised of WT, TBXT heterozygous (TBXT-Het), or TBXT null (TBXT-KO) hiPSCs reveal that varying TBXT dosage does not compromise a cell's ability to differentiate into nascent mesoderm, but that the loss of TBXT significantly delays the temporal progression of the epithelial to mesenchymal transition (EMT). This delay is dependent on TBXT dose, as cells heterozygous for TBXT proceed with EMT at an intermediate pace relative to WT or TBXT-KO. By differentiating iPSCs of the allelic series into nascent mesoderm in a monolayer format, we further illustrate that TBXT dose directly impacts the persistence of junctional proteins and cell-cell adhesions. These results demonstrate that EMT progression can be decoupled from the acquisition of mesodermal identity in the early gastrula and shed light on the mechanisms underlying human embryogenesis.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project.

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

A Mesp1-dependent developmental breakpoint in transcriptional and epigenomic specification of early cardiac precursors.

Transcriptional networks governing cardiac precursor cell (CPC) specification are incompletely understood owing, in part, to limitations in distinguishing CPCs from non-cardiac mesoderm in early gastrulation. We leveraged detection of early cardiac lineage transgenes within a granular single-cell transcriptomic time course of mouse embryos to identify emerging CPCs and describe their transcriptional profiles. Mesp1, a transiently expressed mesodermal transcription factor, is canonically described as an early regulator of cardiac specification. However, we observed perdurance of CPC transgene-expressing cells in Mesp1 mutants, albeit mislocalized, prompting us to investigate the scope of the role of Mesp1 in CPC emergence and differentiation. Mesp1 mutant CPCs failed to robustly activate markers of cardiomyocyte maturity and crucial cardiac transcription factors, yet they exhibited transcriptional profiles resembling cardiac mesoderm progressing towards cardiomyocyte fates. Single-cell chromatin accessibility analysis defined a Mesp1-dependent developmental breakpoint in cardiac lineage progression at a shift from mesendoderm transcriptional networks to those necessary for cardiac patterning and morphogenesis. These results reveal Mesp1-independent aspects of early CPC specification and underscore a Mesp1-dependent regulatory landscape required for progression through cardiogenesis.

Graded mesoderm assembly governs cell fate and morphogenesis of the early mammalian heart.

Using four-dimensional whole-embryo light sheet imaging with improved and accessible computational tools, we longitudinally reconstruct early murine cardiac development at single-cell resolution. Nascent mesoderm progenitors form opposing density and motility gradients, converting the temporal birth sequence of gastrulation into a spatial anterolateral-to-posteromedial arrangement. Migrating precardiac mesoderm does not strictly preserve cellular neighbor relationships, and spatial patterns only become solidified as the cardiac crescent emerges. Progenitors undergo a mesenchymal-to-epithelial transition, with a first heart field (FHF) ridge apposing a motile juxta-cardiac field (JCF). Anchored along the ridge, the FHF epithelium rotates the JCF forward to form the initial heart tube, along with push-pull morphodynamics of the second heart field. In Mesp1 mutants that fail to make a cardiac crescent, mesoderm remains highly motile but directionally incoherent, resulting in density gradient inversion. Our practicable live embryo imaging approach defines spatial origins and behaviors of cardiac progenitors and identifies their unanticipated morphological transitions.

Smarcd3 is an epigenetic modulator of the metabolic landscape in pancreatic ductal adenocarcinoma.

Pancreatic cancer is characterized by extensive resistance to conventional therapies, making clinical management a challenge. Here we map the epigenetic dependencies of cancer stem cells, cells that preferentially evade therapy and drive progression, and identify SWI/SNF complex member SMARCD3 as a regulator of pancreatic cancer cells. Although SWI/SNF subunits often act as tumor suppressors, we show that SMARCD3 is amplified in cancer, enriched in pancreatic cancer stem cells and upregulated in the human disease. Diverse genetic mouse models of pancreatic cancer and stage-specific Smarcd3 deletion reveal that Smarcd3 loss preferentially impacts established tumors, improving survival especially in context of chemotherapy. Mechanistically, SMARCD3 acts with FOXA1 to control lipid and fatty acid metabolism, programs associated with therapy resistance and poor prognosis in cancer. These data identify SMARCD3 as an epigenetic modulator responsible for establishing the metabolic landscape in aggressive pancreatic cancer cells and a potential target for new therapies.

Cell Layers: uncovering clustering structure in unsupervised single-cell transcriptomic analysis.

Motivation: Unsupervised clustering of single-cell transcriptomics is a powerful method for identifying cell populations. Static visualization techniques for single-cell clustering only display results for a single resolution parameter. Analysts will often evaluate more than one resolution parameter but then only report one. Results: We developed Cell Layers, an interactive Sankey tool for the quantitative investigation of gene expression, co-expression, biological processes and cluster integrity across clustering resolutions. Cell Layers enhances the interpretability of single-cell clustering by linking molecular data and cluster evaluation metrics, providing novel insight into cell populations. Availability and implementation: https://github.com/apblair/CellLayers.

Transcription Factor GATA4 Regulates Cell Type-Specific Splicing Through Direct Interaction With RNA in Human Induced Pluripotent Stem Cell-Derived Cardiac Progenitors.

BACKGROUND: GATA4 (GATA-binding protein 4), a zinc finger-containing, DNA-binding transcription factor, is essential for normal cardiac development and homeostasis in mice and humans, and mutations in this gene have been reported in human heart defects. Defects in alternative splicing are associated with many heart diseases, yet relatively little is known about how cell type- or cell state-specific alternative splicing is achieved in the heart. Here, we show that GATA4 regulates cell type-specific splicing through direct interaction with RNA and the spliceosome in human induced pluripotent stem cell-derived cardiac progenitors. METHODS: We leveraged a combination of unbiased approaches including affinity purification of GATA4 and mass spectrometry, enhanced cross-linking with immunoprecipitation, electrophoretic mobility shift assays, in vitro splicing assays, and unbiased transcriptomic analysis to uncover GATA4's novel function as a splicing regulator in human induced pluripotent stem cell-derived cardiac progenitors. RESULTS: We found that GATA4 interacts with many members of the spliceosome complex in human induced pluripotent stem cell-derived cardiac progenitors. Enhanced cross-linking with immunoprecipitation demonstrated that GATA4 also directly binds to a large number of mRNAs through defined RNA motifs in a sequence-specific manner. In vitro splicing assays indicated that GATA4 regulates alternative splicing through direct RNA binding, resulting in functionally distinct protein products. Correspondingly, knockdown of GATA4 in human induced pluripotent stem cell-derived cardiac progenitors resulted in differential alternative splicing of genes involved in cytoskeleton organization and calcium ion import, with functional consequences associated with the protein isoforms. CONCLUSIONS: This study shows that in addition to its well described transcriptional function, GATA4 interacts with members of the spliceosome complex and regulates cell type-specific alternative splicing via sequence-specific interactions with RNA. Several genes that have splicing regulated by GATA4 have functional consequences and many are associated with dilated cardiomyopathy, suggesting a novel role for GATA4 in achieving the necessary cardiac proteome in normal and stress-responsive conditions.

Modeling congenital heart disease: lessons from mice, hPSC-based models, and organoids.

Congenital heart defects (CHDs) are among the most common birth defects, but their etiology has long been mysterious. In recent decades, the development of a variety of experimental models has led to a greater understanding of the molecular basis of CHDs. In this review, we contrast mouse models of CHD, which maintain the anatomical arrangement of the heart, and human cellular models of CHD, which are more likely to capture human-specific biology but lack anatomical structure. We also discuss the recent development of cardiac organoids, which are a promising step toward more anatomically informative human models of CHD.

Targeting transcription in heart failure via CDK7/12/13 inhibition.

Heart failure with reduced ejection fraction (HFrEF) is associated with high mortality, highlighting an urgent need for new therapeutic strategies. As stress-activated cardiac signaling cascades converge on the nucleus to drive maladaptive gene programs, interdicting pathological transcription is a conceptually attractive approach for HFrEF therapy. Here, we demonstrate that CDK7/12/13 are critical regulators of transcription activation in the heart that can be pharmacologically inhibited to improve HFrEF. CDK7/12/13 inhibition using the first-in-class inhibitor THZ1 or RNAi blocks stress-induced transcription and pathologic hypertrophy in cultured rodent cardiomyocytes. THZ1 potently attenuates adverse cardiac remodeling and HFrEF pathogenesis in mice and blocks cardinal features of disease in human iPSC-derived cardiomyocytes. THZ1 suppresses Pol II enrichment at stress-transactivated cardiac genes and inhibits a specific pathologic gene program in the failing mouse heart. These data identify CDK7/12/13 as druggable regulators of cardiac gene transactivation during disease-related stress, suggesting that HFrEF features a critical dependency on transcription that can be therapeutically exploited.

Transcription factor protein interactomes reveal genetic determinants in heart disease.

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

Brahma safeguards canalization of cardiac mesoderm differentiation.

Differentiation proceeds along a continuum of increasingly fate-restricted intermediates, referred to as canalization1,2. Canalization is essential for stabilizing cell fate, but the mechanisms that underlie robust canalization are unclear. Here we show that the BRG1/BRM-associated factor (BAF) chromatin-remodelling complex ATPase gene Brm safeguards cell identity during directed cardiogenesis of mouse embryonic stem cells. Despite the establishment of a well-differentiated precardiac mesoderm, Brm-/- cells predominantly became neural precursors, violating germ layer assignment. Trajectory inference showed a sudden acquisition of a non-mesodermal identity in Brm-/- cells. Mechanistically, the loss of Brm prevented de novo accessibility of primed cardiac enhancers while increasing the expression of neurogenic factor POU3F1, preventing the binding of the neural suppressor REST and shifting the composition of BRG1 complexes. The identity switch caused by the Brm mutation was overcome by increasing BMP4 levels during mesoderm induction. Mathematical modelling supports these observations and demonstrates that Brm deletion affects cell fate trajectory by modifying saddle-node bifurcations2. In the mouse embryo, Brm deletion exacerbated mesoderm-deleted Brg1-mutant phenotypes, severely compromising cardiogenesis, and reveals an in vivo role for Brm. Our results show that Brm is a compensable safeguard of the fidelity of mesoderm chromatin states, and support a model in which developmental canalization is not a rigid irreversible path, but a highly plastic trajectory.

Co-emergence of cardiac and gut tissues promotes cardiomyocyte maturation within human iPSC-derived organoids.

During embryogenesis, paracrine signaling between tissues in close proximity contributes to the determination of their respective cell fate(s) and development into functional organs. Organoids are in vitro models that mimic organ formation and cellular heterogeneity, but lack the paracrine input of surrounding tissues. Here, we describe a human multilineage iPSC-derived organoid that recapitulates cooperative cardiac and gut development and maturation, with extensive cellular and structural complexity in both tissues. We demonstrate that the presence of endoderm tissue (gut/intestine) in the organoids contributed to the development of cardiac tissue features characteristic of stages after heart tube formation, including cardiomyocyte expansion, compartmentalization, enrichment of atrial/nodal cells, myocardial compaction, and fetal-like functional maturation. Overall, this study demonstrates the ability to generate and mature cooperative tissues originating from different germ lineages within a single organoid model, an advance that will further the examination of multi-tissue interactions during development, physiological maturation, and disease.

Dissecting CTCF site function in a tense HoxD locus.

In this issue of Genes & Development, Amândio and colleagues (pp. 1490-1509) dissect the function of a cluster of several CTCF binding sites in the HoxD cluster by iterative deletions in mice. They found additive functions for some, and intriguingly found that some sites function as insulators, while others function as anchors for enhancer-promoter interactions. These functions vary depending on developmental context. The work provides new insights into the roles played by CTCF in regulating developmental patterns and 3D chromatin organization.

Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency.

[Figure: see text].

Modeling Human TBX5 Haploinsufficiency Predicts Regulatory Networks for Congenital Heart Disease.

Haploinsufficiency of transcriptional regulators causes human congenital heart disease (CHD); however, the underlying CHD gene regulatory network (GRN) imbalances are unknown. Here, we define transcriptional consequences of reduced dosage of the CHD transcription factor, TBX5, in individual cells during cardiomyocyte differentiation from human induced pluripotent stem cells (iPSCs). We discovered highly sensitive dysregulation of TBX5-dependent pathways-including lineage decisions and genes associated with heart development, cardiomyocyte function, and CHD genetics-in discrete subpopulations of cardiomyocytes. Spatial transcriptomic mapping revealed chamber-restricted expression for many TBX5-sensitive transcripts. GRN analysis indicated that cardiac network stability, including vulnerable CHD-linked nodes, is sensitive to TBX5 dosage. A GRN-predicted genetic interaction between Tbx5 and Mef2c, manifesting as ventricular septation defects, was validated in mice. These results demonstrate exquisite and diverse sensitivity to TBX5 dosage in heterogeneous subsets of iPSC-derived cardiomyocytes and predicts candidate GRNs for human CHDs, with implications for quantitative transcriptional regulation in disease.

WAPL maintains a cohesin loading cycle to preserve cell-type-specific distal gene regulation.

The cohesin complex has an essential role in maintaining genome organization. However, its role in gene regulation remains largely unresolved. Here we report that the cohesin release factor WAPL creates a pool of free cohesin, in a process known as cohesin turnover, which reloads it to cell-type-specific binding sites. Paradoxically, stabilization of cohesin binding, following WAPL ablation, results in depletion of cohesin from these cell-type-specific regions, loss of gene expression and differentiation. Chromosome conformation capture experiments show that cohesin turnover is important for maintaining promoter-enhancer loops. Binding of cohesin to cell-type-specific sites is dependent on the pioneer transcription factors OCT4 (POU5F1) and SOX2, but not NANOG. We show the importance of cohesin turnover in controlling transcription and propose that a cycle of cohesin loading and off-loading, instead of static cohesin binding, mediates promoter and enhancer interactions critical for gene regulation.